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We present a systematic investigation of molecule-metal interactions for transition-metal phthalo- 
cyanines (TMPc, with TM = Fe, Co, Ni, Cu) adsorbed on Ag(lOO). Scanning tunneling spectroscopy 
and density functional theory provide insight into the charge transfer and hybridization mechanisms 
of TMPc as a function of increasing occupancy of the 3d metal states. We show that all four TMPc 
receive approximately one electron from the substrate. Charge transfer occurs from the substrate 
to the molecules, inducing a charge reorganization in FePc and CoPc, while adding one electron 
to ligand 7r-orbitals in NiPc and CuPc. This has opposite consequences on the molecular magnetic 
moment: in FePc and CoPc the interaction with the substrate tends to reduce the TM spin, whereas 
in NiPc and CuPc an additional spin is induced on the aromatic Pc ligand, leaving the TM spin un- 
perturbed. In CuPc, the presence of both TM and ligand spins leads to a triplet ground state arising 
from intramolecular exchange coupling between d and tt electrons. In FePc and CoPc the magnetic 
moment of C and N atoms is antiparallel to that of the TM. The different character and symmetry 
of the frontier orbitals in the TMPc series leads to varying degrees of hybridization and correlation 
effects, ranging from the mixed- valence (FePc, CoPc) to the Kondo regime (NiPc, CuPc). Coherent 
coupling between Kondo and inelastic excitations induces finite-bias Kondo resonances involving 
vibrational transitions in both NiPc and CuPc and triplet-singlet transitions in CuPc. 
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I. INTRODUCTION 

The interaction of molecules with solid s urfac es governs 
the self-assembly of supramolecular layerJilll as well as 
the electronic and magnetic properties of metal -org anic 
heterostructurea^El and single molecule devices.'^ The 
large variety of organic and metal-organic complexes, 
combined with the possibility of tuning the chemical reac- 
tivity and electronic ground state of different complexes 
within the same family of molecules, make such systems 
increasingly attractive for applications. However, elec- 
tronic interactions at the interface with a metal may in- 
duce charge transfer, distort the ligand field, and reduce 
electron-electron correlation effects via screening and hy- 
bridization. The intricate interplay between all these pro- 
cesses is far from being understood. 

Transition-metal phthalocyanines (TMPc) represent a 
well-known class of molecules with organic semiconduct- 
ing properties in the bulk and a broad spectrum of ap- 
plications that includes field-effect-transistors, gas sen- 
sors, and photovoltaic cells.'^ Due to their relatively sim- 
ple and robust structure as well as versatile chemistry, 
TMPc have assumed the role of a model system to study 
the interaction of metal-organic complexes with metal 
surfaces. Their flat adsorption geometry facilitates the 
bonding of both the central TM ion and organic ligands 
to the substrate, whereas their capability to coordinate 
many different metal atoms allows for a systematic in- 



vestigation of their magnetic properties. 

The interaction of TMPc with metal surfaces has been 
studied by many techniques, including scanning tun- 
neling microscopy (STM),'2EH25] ultrav iolet and X-ray 
photoelectron spectroscopy ( UPS-X PS) jMUMl x-ray ab- 
sorption spectroscopy (X AS), I^EIMtll] and density func- 
tional theory (DFT).'22H2ni investigations of the ma gnetic 
prop erties of T MPc have mostly fo cused on MnPcP^I^, 
FePclSIHESMl and CoPcl^ElHiailHSl] ^j. ^^ich it was 

shown that adsorption on a metallic surface tends to 
quench the magnetic moment of the TM ion. A corol- 
lary to such investigations is the implicit assumption 
that the magnetic properties of TMPc at surfaces de- 
pend almost exclusively on the ground state of the TM 
ion. However, a recent investigation by our group has 
highlighted the possibility of inducing a magnetic mo- 
ment delocalized over the Pc organic ligand, which af- 
fects the spin degeneracy and symmetry of the molec- 
ular ground state as well as the electrical conductance 
measured at different sites within the same molecule.'^ 
Unpaired spins induced by charge transfer in Pc ligands 
have also been found to affect the conductance proper- 
ties of TbPc2 single molecule magnets adsorbed on AiP^l 
and their coupling to ferromagnetic surfaces.!^ Indeed, 
a growing number of experiments indicates that the or- 
ganic ligand directly affects th e magn etism and trans- 
port properties of metal-organicpSESEfll as well as purely 
organic complexe^^J^idsorbed on surfaces. 

These results evidence the need to take into account 



the full molecular structure in order to understand and 
predict the magnetic properties of adsorbed systems. 
Here we present a comprehensive study of four differ- 
ent TMPC (TM = Fe, Co, Ni, Cu) deposited on Ag(lOO) 
in order to shed light on the mechanisms that lead to 
changes of both the TM and ligand magnetic moments 
at the TMPc/metal interface. Scanning tunneling spec- 
troscopy (STS) combined with DFT allows us to eluci- 
date the role of charge transfer, hybridization, and cor- 
relation at the molecule-metal interface on the magnetic 
properties. The paper is organized as follows: after pre- 
senting the experimental and theoretical methods (Sec- 
tionln]) , we first introduce the theoretical electronic struc- 
ture of TMPc in the gas-phase (Section |lII[ ), which serves 
as the starting point of our investigation. We then de- 
scribe the experimental and theoretical adsorption geom- 
etry of the molecules on Ag(lOO) (Section IV), including 



a discussion of electronic chirality effects, which were re- 
vealed in a previous work on CuPc.^ The electronic struc- 
ture of adsorbed TMPc is discussed in Section [V| where 
STS data are compared with the molecular projected 
density of states (PDOS) derived from ab-initio DFT 
calculations. The magnetic properties of TMPc are ad- 



dressed in Section VI using STS to measure the intensity 
and spatial extension of Kondo anomalies in the differ- 
ential conductance, including the coupling to vibrational 
and magnetic degrees of freedom in NiPc and CuPc as 
well as the theoretical description of triplet-singlet exci- 
tations in CuPc using a multi-orbital approach. Further, 
we calculate the magnetic moments of adsorbed TMPc 
and compare them to the gas-phase magnetic moments 
in neutral and negatively charged species in order to un- 
derstand the effects of charge transfer and intramolecular 
exchange and their dependence on the TM ion. Finally, 
Section [VII| summarizes the more general conclusions of 
this work. 



II. METHODS 



A. Experiment 



tip or substrate electronic structure.ESl Hence, only rela- 
tive dl/dV values are significant in these spectra. dl/dV 
maps were acquired in the constant current mode, with 
an amplitude of 3 mVrms ■ Maps of the Kondo resonance 
were obtained using the d^I/dV^ signal in order to scan 
at a finite voltage-l^ Both topographic and conductance 
images were processed using the WSxM software. Fur- 
ther details on the spectroscopic methods can be found 
in the Supplementary Information of ReflT] 



B. Theory 

The theoretical electronic structure of both gas- 
phase and adsorbed TMPc have been obtained from 
ab-initio calculations using the VASP implementa- 
tion of DFT in the projected augmented plane wave 
scheme.'^^ESl Differe nt ap proximations, namely LDA,SSl 
GGA,ElGGA-hvdW,l3SllS!and GGA+lPlhave been com- 
pared in order to obtain consistent results, and to check 
the effect of different exchange-correlation approxima- 
tions or the vdW interaction. Unless otherwise stated, 
the results presented in the next sections were obtained 
using GGA+l^ and vdW corrections using the scheme 
of Grimme.SSBSl Ugj^=U-J was chosen to be 3 eV in all 
cases. Additional tests were made to check that reason- 
able values of Ue// do not qualitative change our conclu- 
sions. The planewave cutoff energy was set to 300 eV. 
The calculated slab included 5 Ag atomic layers interca- 
lated by 7 vacuum layers in the vertical direction, and a 
7x7 lateral supercell. The positions of all atoms in the 
molecule and the first three Ag layers were relaxed ver- 
tically and laterally until forces were smaller than 0.05 
eV/A. The projected density of states (PDOS) of NiPc 
have been used to compare single (1x1x1) and multi- 
ple (5x5x1) k point calculations, which converge after 
applying a broadening of 100 meV to the data. Based 
on that, the results for a single k point with the latter 
broadening are used in the following. Charge transfer 
and local magnetic moments have been calculated using 
a Bader charge analysis. ^-1^^ 



TMPc were evaporated in ultra high vacuum (UHV) 
from a heated quartz crucible onto a sputter-annealed 
Ag(lOO) single crystal kept at room temperature, after 
degassing the 99% pure powder material (Sigma Aldrich) 
for 24 hours. The deposition rate was ~ 0.05 monolay- 
ers/min; the base pressure during evaporation was below 
5 X 10"^" mbar. Spectroscopic measurements were per- 
formed using an STM operating at 5 K. Conductance 
{dl/dV) spectra were obtained with the lock-in tech- 
nique, using a bias voltage modulation of frequency 3 
kHz and amplitude 1 mVrms for the low energy spectra 
(—0.1 < V < 0.1 V), and 3 mV^ms for the larger energy 
range (—2 <V<1 V). A background spectrum acquired 
on the bare Ag surface with the same feedback conditions 
was subtracted to all spectra in order to enhance molecu- 
lar resonances and minimize features originating from the 



III. GAS-PHASE TMPC 

Before discussing the effects of adsorption on the 
molecular electronic and magnetic structure, we briefly 
introduce the chemical configuration of gas-phase TMPc, 
as well as the electronic structure obtained from spin- 
polarized calculations, which will be used later as a ref- 
erence. 

TMPc are characterized by a rather simple and sta- 
ble molecular structure, where the TM ion bonds to four 
isoindole (pyrrole -I- benzene) ligands via their pyrrole N 
(Np). That leaves the ion in a [TM]^+ state, and the 
molecule with the square planar D^h symmetry. Under 
this symmetry group, the d-states transform as 623 (d^y), 
big {d^2_y2), aig (42), and Cg {dxz4yz = d-^). Depend- 
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FIG. 1; (Color online). Spin-polarized electronic structure of gas-phase TMPc. (Un-) occupied states are represented by 
(dashed) full lines. The contribution of different d-states is indicated in percentage. 



ing on their symmetry and energy position, these orbitals 
mix to a different degree with 2p-states of the C and 
N atoms. The highest occupied and lowest unoccupied 
molecular orbitals (HOMO and LUMO) of the Pc ring are 
represented by delocalized ai„ and 2eg 7r-orbitals, respec- 
tively, with marginal contribution from the TM d-states. 
The TM-related MO can be classified in two groups ac- 
cording to their parallel (&2g, big = dy) or perpendicular 
(aig, Cg — d±) orientation with respect to the molecular 
plane, and hence, the substrate. 

Figure [T] shows the spin-up and spin-down energy lev- 
els of TMPc near the frontier orbitals obtained within 
the GGA+U scheme. The 3d contribution in each MO 
is indicated between parentheses; the molecular magnetic 
moments are discussed in Sect. IVfl Our calculations show 
that the electronic structure of the Pc ring is barely af- 
fected by the TM ion. This is clearly seen in the aiu 
and 2eg ligand orbitals: the aiu MO, with negligible d- 
contribution, appears at the same energy in all molecules. 
The 2eg energy varies slightly with the d configuration of 
the TM ion as this orbital exhibits a finite d^ contribu- 
tion that ranges from 7% in FePc to 1% in CuPc. 

The evolution of the electronic configuration of the cen- 
tral [TM]^+ ions cannot be rationalized by the simple 
filling of one-electron energy levels split by the ligand 
field, due to the strong d-d correlation and the energy- 
dependent hybridization with ligand states. The sit- 
uation is most critical for the ground state of FePc, 
where the quasidegeneracy between different electronic 



terma^ESl has led to a long-standing debate in litera- 
|-yj.gj54H59] rj^Yie discussion of such differences is beyond 
the scope of this paper. However, the spin moment that 
we obtain for this molecule is 8=1, independently of the 
DFT fun ctional and in good agreement with previous 
work.'SEMSSl CoPc, with one electron more, presents a ro- 
bust ^Aig (S=l/2) ground state with a single aig (^22) 
hole, in agreement with x-ray absorption measurements.'^l 
NiPc, with d^, has a singlet (S=0) ^Aig ground state. 
The TM spin is recovered in CuPc, with a single hole in 
big {dx2_y2), resulting in a doublet (S=l/2) ^Big ground 
state. 



IV. ADSORPTION OF TMPC ON AG(IOO) 

A. Adsorption geometry 

The adsorption of individual TMPc has been stud- 
ied by STM. Their topographic appearance already re- 
flects differences that allow us to classify them in two 
groups: the TM ion appears as a protrusion in FePc and 
CoPc and a depression in NiPc and CuPc, as shown in 
Figs. [2] (a) and [SJ The different contrast is due to the 
large (small) coupling of the TM d±_ {cIm) orbitals near 
the Fermi level {Ep) with the tip states.!^ As typical for 
TMPc on metals, the molecules adsorb with the aromatic 
plane parallel to the surface. Atomically resolved images 
of the surface reveal two different azimuthal orientations 



of the TM-Np ligand axis [Fig. [2] (b)], corresponding 
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to ±30° with respect to the (lOOj surface lattice vec- 
tors. Previous systematic ab-initio calculations of CuPc 
on Ag(lOO) revealed that this orientation, with the TM 
ion located on the Ag hollow site, corresponds to the min- 
imum of the adsorption energy-l^l A similar adsorption site 
has been found for CuPc on Cu(lOO) (Ref.EI]). The ro- 
tation of the Pc macrocycle with respect to the substrate 
high symmetry directions was attributed to bond opti- 
mization between aza-N (Na) and Ag atoms. This study 
confirms the leading role of the N^-Ag interaction, as 
both theory and experiment find very similar adsorption 
configurations for all TMPc, independently of the central 
TM ion (see Fig. [2})). 

The calculations also show that TMPc become slightly 
concave after deposition, with the periphery of the 
molecule slightly further from the surface. This effect 
is most pronounced in the case of FePc. The molecules 
also induce a small distortion of the substrate, push- 
ing the Ag atoms below Np slightly below the surface 
plane and the Ag atoms below Na slightly above it. The 
molecule-substrate distances, calculated as the difference 
in the z coordinate between the TM atom and the Ag 
atom below Np, are shown in Table ll| The distance ob- 
tained with LDA is ^-^2.5 A for the four molecules. This 
method is known to overbind, compensating for the ab- 
sence of van der Waals (vdW) forces. These forces are 
included in the GGA-|-vdW method. The results with 
GGA+vdW show that z increases to ~2.7 A, keeping it 
TM-independent. Similar values were also found using 
both methods for CoPc adsorbed on Gu(lll) (Ref. 134)1 . 
On the less reactive Au(lll) surface, however, distances 
obtained with LDA for different TMPc do not le vel out 
and still reflect a TM-dependent behavior.lSSEn! jn our 
case, we can separate the effect of the interaction between 
the TM ion and the substrate by using plain GGA, i.e., 
by switching off the vdW interaction. Using this method, 
we find that FePc and CoPc, with z = 2.76 and 3.08 A 
respectively, are significantly closer to the surface than 
NiPc and CuPc, with z = 3.59 and 3.47 A respectively. 
The stronger TM-Ag interaction of FePc and CoPc is 
attributed to the direct partici pation o f the d-orbitals 
to the molecule-substrate bonds!2£l22l36| The vdW inter- 
action introduced by GGA-|-vdW can be estimated by 
the correction added to the total energy calculated by 
GGA alone. This correction is larger by about 1 eV for 
NiPc and CuPc compared to FePc, which compensates 
for the lack of direct TM-substrate interaction in the for- 
mer species and results in the nearly TM-independent 
bond length reported in Table |l) 



B. Electronic chirality 

Figure [3] shows that NiPc and CuPc present chiral con- 
trast at negative bias, whereas FePc and CoPc appear 
achiral at both negative and positive bias. This effect has 
been studied in detail for CuPc/Ag(100) in Ref. 1^, where 
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FIG. 2: (Color online), a) Topographic image of CoPc 
and NiPc molecules co-deposited at room temperature on 
the Ag(lOO) surface (7 = 0.49 nA, V = -1.0 V), image size 
13.4 X 13.4 nm^. The green/red arrow indicates the +/ — 30° 
rotation of the molecular axis with respect to the surface lat- 
tice vectors (white arrows), b) Adsorption geometry of NiPc 
on Ag(lOO) calculated by DFT. The different elements are la- 
beled. Ni and N2 refers to aza and pyrrole N, and Ca and 
Cs to the two C atoms of the benzene ring used to measure 
chiral distortions (Table ^. 
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TABLE L Computed height distances z (in A) between the 
TM ion and the Ag atom below Np, obtained within different 
DFT approximations, and height difference between Ca and 
Cs atoms at the end of the benzene ring (see Fig.[2p) obtained 
with GGA-hvdW. 



it has been shown that the chiral contrast observed by 
STM is not related to nuclear distortions, but originates 
from the asymmetric electronic interaction between the 
aiu and 2eg orbitals and Ag states. This asymmetry is 
due to the misalignment between surface and molecular 
symmetry axes, where the -1-30° and —30° configurations 
imprint chiral charge distortion of opposite sign. Here we 
extend this study to FePc, CoPc, and NiPc, and quan- 
tify the contribution of conformal distortions to the chi- 
ral contrast. We use the height difference between two 
opposite C atoms at the benzene ring as a measure of 
their torsion (see Table |I] and Fig. [2]d). Its small value, 
close to the calculation's accuracy, obtained for all TMPc 
confirms the electronic origin of chirality. This is in line 
with the voltage dependence of the chiral contrast, and 
the fact that it is only observed in CuPc and NiPc in 
spite of the similar geometry found for all molecules. 

Electronic chirality is imprinted to each MO to a dif- 
ferent degree, depending on their spatial distribution and 




FIG. 3: (Color online). Bias dependent topographic images 
of all TMPc. For CuPc and NiPc negative voltage images 
with maximum chiral contrast are displayed. The contrast of 
CoPc and FePc does not vary in the range of negative voltages 
investigated. For positive voltage, all TMPc appear achiral 
in the bias range V < +1 V. 



hybridization with the substrate, with maximuni effect 
in orbitals exhibiting a nodal plane on the TM-Np ligand 
axis (see Fig. |5|. The ai„ MO exhibits the strongest 
distortion because of its double lobe structure at the 
benzene ring, where the asymmetry of charge distortion 
is maximum. The chiral contrast at the 2eg MO, also 
with nodal planes on the ligand axis, varies with en- 
ergy: at positive bias it is strongly reduced, but becomes 
more intense when shifted to negative bias due to charge 
transfer. The MOs with significant TM-d contribution 
do not have nodal planes on the ligand axis, and hence 
do not present any significant chirality. This leads to 
a voltage-dependent chiral appearance of the TMPc that 
depends solely on the specific electronic structure of each 
molecule, as shown in Fig.|3] CuPc and NiPc are strongly 
chiral at negative bias due to the asymmetry of the ai„ 
orbital and of the partially occupied 2eg orbital, and achi- 
ral at positive bias, where tunneling occurs via a more 
symmetric 2eg resonance. In CoPc and FePc, on the 
other hand, the 2eg orbital only appears at positive bias, 
where the chiral contrast is insignificant. Furthermore, 
TM-d states and other other interface-related resonances 
with maximum intensity on the ligand axis hinder the 
chirality of the ai„ state. 



V. ELECTRONIC STRUCTURE 

In this section we present scanning tunneling spec- 
troscopy (STS) measurements of the molecules adsorbed 
on the Ag(lOO) surface. The energy of the MOs and 
other interface-related features is obtained from individ- 
ual dl/dV spectra, whereas maps of the dl/dV intensity 
at specific energies reveals their symmetry and spatial 
distribution within the molecule. Experimental results 
are then compared to the theoretical density of states 
projected onto different MOs, which allows us to assign 
the observed dl/dV peaks to specific orbitals, provid- 
ing insight into the formation of hybridized molecule- 
substrate states, charge transfer, and molecular mag- 



netism. 



A. Scanning tunneling spectroscopy 

Figure Qa) shows the dl/dV spectra of FePc to CuPc 
recorded by positioning the STM tip on the central TM 
atom (red curves) and on a peripheral benzene ring (blue 
curves). The aiu and 2eg resonances corresponding to 
the HOMO and LUMO of the gas-phase Pc ring can be 
easily identified by studying the spatial distribution of 
the dl /dV peaks of the spectra acquired on benzene, as 
shown in Fig. [sj The intensity distribution of these res- 
onances is similar to that found for other substrates,'2Sl 
except for the chiral fingerprint discussed in the previous 
Section. The energy position of the aiu state is simi- 
lar in all TMPc, varying continuously from -1.14 V in 
CoPc to —1.40 V in CuPc (FePc was not measured in 
this energy range). On the other hand, the spectral dis- 
tribution of the 2eg state critically depends on the TM 
ion: in FePc and CoPc, a single unoccupied peak is ob- 
served at -1-0.39 and -1-0.47 V, respectively, as expected 
for the Pc LUMO. Conversely, CuPc and NiPc exhibit 
two peaks around the Fermi level, located at -0.29/-f 0.35 
V and -0.35/-I-0.35 V, respectively. Both peaks, however, 
can be associated with the 2eg state due to their similar 
spatial distribution. The position of the lowest energy 
peak below Ep indicates that this orbital is partially oc- 
cupied. The energy splitting of ~ 0.65 — 0.70 eV is within 
the range of the Coulomb repulsion energies obtai ned for 
TT-orbitals in aromatic complexes of similar size.'HMM] 
We therefore assign the peaks at negative and positive 
bias to single and double occupation of the 2eg state, 
respectively, in agreement with the results reported in 
Ref. T'. The sharp peak appearing between the two 2eg 
"Coulomb blockade" peaks is the Kondo resonance as- 
sociated to the unpaired ligand spin,'^ which will be dis- 
cussed further in Section IVll Note that the short electron 
lifetime of empty orbitals in molecules coupled to metals 
prevents their double occupation by tunneling electrons, 
which would lead to similar "Coulomb blockade" peaks 
in the 2eg state of FePc and CoPc. 

Spectra recorded on the TM ions, shown as red curves 
in Fig. |4^a), allow us to identify the MOs with d- 
character. The broad peak shifting downwards with in- 
creasing d-occupation from FePc to NiPc can be related 
to (i_L states, easiest to observe due to their effective cou- 
pling with the STM tip. In CuPc, these states are situ- 
ated below the probed energy range. The localization of 
this peak on the TM ion is confirmed by the dl / dV map 
of CoPc at -510 mV [Fig.[5];b)]. In this molecule, the dj_ 
resonance appears well below Ep, suggesting the com- 
plete filling of the aig (^^2) orbital, and hence a switch 
in the charge transfer channel from the ligand 2eg to TM- 
d states compared to NiPc and CuPc. In FePc, the tail of 
the d_L peak crosses Ep, suggesting that these states are 
not fully occupied, similarly to FePc on Au(lll) (Ref.|3|. 

Contrary to the dj^-states, the perturbation of the du- 
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FIG. 4: (Color online), (a) dl/dV spectra, acquired on the TM ions (red, bottom) and benzene rings (blue, top). The latter 
spectra have been offset for clarity. MOs, interface (IR) and Kondo (K) resonances, as well as Fano-type (F) features are 
labeled. Tunneling conditions (/, V): 1 nA, -1.0 V (FePc, CoPc), 3 uA, -2.2 V (NiPc), 3 nA, -2.0 V (CuPc). (b) Computed 
spin-polarized DOS of TMPc (black), and its projection onto the TM d-states. (c) Computed spin-polarized PDOS projected 
onto MOs. 



states due to adsorption appears to be minor. Although 
confinement in the molecular plane makes them difficult 
to access by STS, some features in the spectra of NiPc 
and CuPc can be tentatively assigned to the big (dx'^-y'^) 
state. In the pristine molecules, this orbital is singly 
occupied in CuPc and unoccupied in NiPc (see Fig. IT]). 
The intensity corresponding to the unoccupied big state 
is masked by a step-like interface resonance feature (IRi ) 
that we observe around +0.2 V in all molecules. How- 
ever, both CuPc and NiPc exhibit a larger dl/dV sig- 
nal in this energy range compared to CoPc and FePc, 
which suggests that there is an additional orbital con- 
tributing to the tunneling current above Ep. The singly 
occupied big state in CuPc can be assigned to the small 
hump observed around -0.4 V, close to where the 2eg 
resonance on the benzene spectra presents a pronounced 



valley. This valley, absent in NiPc at this energy, is ten- 
tatively assigned to a Fano resonance originating from 
the tunneling interference between the occupied big and 
the more hybridized 2eg states, analogous to that occur- 
ring between resonant and direct channels in impurity 
and quantum dot systems .'^^HMI Indeed, the dip observed 
in NiPc at -1.4 V, where d± states overlap with the aiu 
orbital, might have a similar origin, supporting the pres- 
ence of Fano peaks in these molecules. 

In addition to MO resonances that maintain the orig- 
inal character of the TMPc states, we find additional in- 
terface resonances that originate from strong hybridiza- 
tion with Ag electrons (see Fig. [6]) . The IRi resonance 
observed at -hO.2 V [Fig. ga)] is common to all TMPc 
and only observed on Ag(lOO). This resonance could orig- 
inate from the hybridization of molecular states with the 




FIG. 5: (Color online). Constant current dl /dV maps of (a) 
CuPc and (b) CoPc. A topographic image of each molecule 
recorded simultaneously during the acquisition of the maps is 
displayed on the left side. The color scale is adjusted for max- 
imum contrast in each map. NiPc and FePc present similar 
MO distributions as CuPc and CoPc, respectively. 




FIG. 6: (Color online). Constant current dl/dV maps of 
interfacial resonances in FePc arising from the hybridization 
between molecular and substrate electrons, marked as IR in 
Fig. |4|b). The topography of FePc recorded during the ac- 
quisition of the maps is displayed on the left side. 



lowest lying unoccupied surface state of Ag(lOO), which 
is close in energy.l^ In FePc, we observe an additional 
feature (IR2) between the ai„ and d± states, which we 
identify with an interface state with intensity at the ben- 
zene rings, similar to that reported for CoPc on Au(lll) 



(Ref.Hni)- 



B. Density of states calculated by DFT 

In order to correlate the features observed in the dl/dV 
spectra with the molecular electronic structure, we have 
computed the PDOS projected onto the TM-d states as 
well as on the MOs of the TMPc, as shown in Fig. |4J In 
general, we find fair agreement between the measured and 
calculated resonances, supporting the assignment of the 
spectral features to specific MO based on the symmetry 
and spatial distribution of the dl/dV intensity. 

The strongest hybridized states are the aig MOs of 
FePc and CoPc, as expected owing to their dominant 
dz2 character. The TM ^22 -electrons hybridize with the 
Ag-spz states, confirming the strong, direct substrate- 



TM interaction discussed in Section |IV] Additionally, in 
the case of FePc and CoPc, two 2eg and one Cg state 
are mixed together, projecting onto one occupied and 
one empty spin-down resonances (see Fig. He). Mixing 
is induced by hybridization with the substrate, showing 
that the picture of slightly distorted MO is no longer 
valid in this case. This is in line with the appearance of 
interface resonance as those observed in STS maps (IR2 
in Fig. l6]). For NiPc and CuPc, hybridization with the 
substrate turns out to be much smaller, and the PDOS 
resembles that of the gas-phase molecules of Fig. [T] In 
these two cases the confinement of the big orbital in the 
molecular plane and its a character reduces hybridization 
with Ag electrons, diminishing the charge transfer to this 
orbital. Hence, it appears unoccupied in NiPc and singly 
occupied in CuPc, as in the gas-phase. 

The aiu orbital, fixed at about -1.30 eV, is the only 
MO that is not affected by the d-occupation of the TM 
ion. In the case of FePc, the spin up and spin down «!„ 
states are exchange-split due to the single-configurational 
nature of DFT. This broken symmetry reflects the strong 
spin-polarization of the whole molecule. The position of 
the 2eg MO, on the other hand, depends strongly on the 
type of TM ion and reveals a transition in the interaction 
with the substrate. As in the dl/dV spectra, the PDOS 
shows that this state is unoccupied in FePc, partially oc- 
cupied in CoPc, and acquires a charge of approximately 
one electron in NiPc and CuPc. The splitting between oc- 
cupied and unoccupied 2eg states is not fully reproduced 
in the calculations due to the well-known electronic gap 
problem in DFT, which is especially critical for dclocal- 
ized TT-orbitals.'SSl 

The origin of the charge-transfer behavior with TM 
can be unfolded by using gas-phase anions to study pure 
transfer without hybridization. We find that the TM- 
d dependent behavior is the same as in adsorbates. In 
[FePc]^~ and [CoPc]^~ the effect of the extra electron is 
complex and the extra charge is not easily assigned to 
a single MO. It rather induces a reorganization of the 
charge within the molecule. In [NiPc]"'^^ and [CuPc]^~ 
the behavior is simpler, with the extra electron going to 
2eg. The observed behavior is also in agreement with 
transport results in electron-doped thin films 
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In summary, although the net amount of transferred 
charge is about one electron in all cases (see Table |ll| , 
its effect on the molecule depends on the d occupation of 
the ion. This difference is driven by the interplay between 
energy alignment and Coulomb repulsion in the d and 
2eg states, as well as their degree of hybridization with 
the substrate. The similarity between results obtained 
by calculating gas-phase anions and adsorbates indicates 
that the leading role of the substrate is that of a charge 
reservoir. However, the strong (weak) hybridization of 
the d± ((i|[) states further stabilizes the different charge 
donation channels, promoting charge transfer to TM-d 
states in FePc and CoPc, and inhibiting it in NiPc and 
CuPc. 



VI. MAGNETIC PROPERTIES 

The magnetic properties of TMPc on Ag(lOO) are ad- 
dressed using STS through the analysis of dl/dV spec- 
tra close to Ep. The presence of sharp zero bias reso- 
nances can be used to detect the formation of local mag- 
netic moments screened by the Kondo effect f^^'^ This 
occurs as itinerant electrons belonging to the substrate 
couple antiferromagnetically to a localized spin, form- 
ing a many-body singlet state at temperatures below the 
characteristic Kondo energy UbTk- dl/dV maps of the 
Kondo spectral features further allow us to measure the 
spatial distribution of the spin density for unpaired elec- 
trons. We note, however, that, although the presence of 
a Kondo resonance is a signature of a local moment, its 
absence is not a proof of a nonmagnetic ground state, 
since the resonance can also be experimentally inacces- 
sible due to the low Tr- in weakly interacting spins. In 
such a case, STS data have to be contrasted with com- 
plementar y ex periments, such as x-ray magnetic circular 
dichroism,SEl! or, as in the present case, with theoretical 
calculations. 
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FIG. 7: (Color online), (a)- (d) Conductance spectra acquired 
around Ep on the TM ions (red) and benzene rings (blue). 
The zero-bias (elastic) , vibrational, and triplet-singlet Kondo 
resonances observed in CuPc and NiPc are labelled as K^, 
Kv, and Kts respectively. Tunneling conditions (7, VJ,): 1 
nA, -100 mV (FePc, CoPc), 1.1 nA, -100 mV (NiPc), 2 nA, 
-100 mV (CuPc). (e) and (f) d^I/dV'^ maps acquired at -3 
mV and -5mV respectively, showing the intensity distribution 
of the elastic Kondo resonance {Ke) of CuPc and NiPc and 
its close resemblance to the 2e„ resonances. 



A. Kondo effect probed by STS 

Figure [t] shows the dl/dV spectra measured near Ep 
for the four molecules at the position of the TM ion 
(red curves) and benzene ligand (blue curves). FePc 
and CoPc present featureless spectra, which indicate that 
the Kondo interaction between these molecules and the 
Ag(lOO) surface is either absent or too weak to be ob- 
served at 5 K. Relatively flat Co spectra were reported 
also for CoPc adsorbed on Au(lll), for which it was con- 
cluded that the filling of the aig state leads to a com- 
plete quenching of the molecular magnetic moment.'^ 
This interpretation has been recently questioned based 
on a mixed- valence model of x-ray absorption spectra.'^ In 
both cases, the quenching of the CoPc magnetic moment 
appears to be a robust result. On the other hand, FePc 
does present a Kondo resonance on Au(lll), alth ough 
different interaction strengths have been reported.'i^E^ 
We attribute the absence of Kondo peaks on Ag(lOO) 
to a stronger interaction with the substrate, as reflected 
by the substantial hybridization of the d± states, which 
may lead to a mixed-valence configuration of FePc, as 
discussed later. 

In contrast with FePc and CoPc, NiPc and CuPc 
present many characteristic signatures of low-energy exci- 
tations around Ep. Such excitations have been discussed 
in detail in Ref. [3 Here we report a summary of these 
results for comparison with theory as well as to provide 
a complete picture of TMPc adsorbed on Ag(lOO). The 
most prominent feature of the NiPc and CuPc spectra 
is the presence of an intense zero bias resonance with 
maximum intensity at the ligand position. This is rather 
unusual for metal-organic adsorbates since, in most cases, 
Kondo resonances are associated to the spin of the TM 



JQjjg jTUTBITglMgs] jjgj-g^ |;jjg 2ero bias peak is assigned to 
an elastic Kondo resonance (Kg) arising from the un- 
paired spin located in the 2eg orbital of CuPc and NiPc. 
This interpretation is supported by the temperature de- 
pendence of the peaks full width at half maximum (Tj^) 
and maximum intensity {Gk)- Figure [s] shows that both 
parameters follow expressions derived from Fermi liquid 
theorjED and numerical renormalization group (NRG) P^ 



VKiT) - 2^{TikBTf + 2{kBTKf 
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Ga-(T) = Go// + G(0) 
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where a and ^ are parameters that depend on the total 
spin of the molecule. Fitting the values of Tk with Eq. [T] 
leads to a Kondo temperature of Tk = 27 ± 2 K for 
CuPc. The logarithmic behavior of Gk makes it more 
sensitive to the scattering of the data, but we find that 
expressions for both S=l/2 and (the underscreened) S=l 
fit reasonably well the data by using Tk — 27 K. The two 
possible spin states will be discussed later. 

We notice that NiPc and CuPc present very similar 
Kondo resonances in spite of the S=0 ground state of the 
Ni ion and S=l/2 ground state of Cu. Indeed, we esti- 
mate a Kondo temperature Tk = 29 K for NiPc from Vk 
measured at 5 K. Moreover, the intensity of both CuPc 
and NiPc resonances closely follows the contour of the 2eg 
orbital, as revealed by the remarkable similarity of the 
dl/dV maps reported in Figs. JTJJe) and (f ) compared to 
Fig. [5] These observations indicate the presence of an un- 
paired spin localized on the 2eg ligand state, in agreement 
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FIG. 8: (Color online). Temperature dependence of (a) Tk 
and (b) Gk normalized to the intensity at 5 K (Go) of the 
elastic Kondo resonance of CuPc. The solid line in (a) is a fit 
of Tk using Eq. [l] giving Tk = 27 ± 2 K. The solid (dashed) 
line in (b) is a fit of Gk /Go for the S=l/2 (underscreened 
S=l) case of Eq. [2] 



B. Inelastic vibrational and Kondo excitations 

The dl/dV spectra of NiPc and CuPc recorded on the 
TM ions [red curves in Figs.JTJJc) and (d)] present a com- 
plex series of peaks and conductance steps in addition to 
the elastic Kondo resonance observed at zero bias. It is 
well-known that inelastic excitations induce step-like in- 
creases in the differential conductance s pectra symmet- 
rically distributed in energy around Ep^^'^ In Kondo 
systems, the coherent coupling of the Kondo state with 
such excitations leads to additional peaks at the ener- 
gies of the inelastic conductance steps.^^ Unlike quantum 
dots, where the reduced orbital level-spacing favors the 
coupling of Kondo and electronic excitations,'^^ the ori- 
gin of such side peaks in sm all aromatic moleculesis re- 
strained to vibrationaPiEsHlo] (^j^^^ gj. magneticPlZaiMiH 

(Kts) excitations. In the present case, by comparing the 
energy of the K^ peaks with that of the Raman vibra- 
tional modes of gas-phase CuPc,'S31 we find that the in- 
elastic features found at the position of the Ni and Cu 
ions [Fig.loja)] correspond to vibrational excitations that 
involve distortion of the TM-Np bonds. ^ Note that the 
vibrational spectra of NiPc and CuPc are nearly iden- 
tical, as shown by the twin cPl/dV'^ spectra displayed 
in the bottom graph. The step-like and peak-like con- 
tributions to the inelastic conductance can be obtained 
by fitting the dl /dV spectra with step and Lorentzian 
functions, as shown in Fig.loFa). From the fit of the spec- 
tra we observe that both NiPc and CuPc present more 
intense conductance steps at negative bias and more in- 
tense Kondo peaks at positive bias. These inverted cor- 
relation between vibrational and Kondo features has also 
been observed for TCNQ molecules adsorbed on Au(lll), 
and has been attributed to the competition between the 
purely inelastic channel and the one including the Kondo 
effect .El 



with the single occupancy of this orbital discussed in Sec- 
tionlV] The relation between energy, width, and Coulomb 
repulsion potential, with values of e = —0.30 eV, F ~ 0.30 
eV, and U = 0.65 eV, respectively, is also characteristic 
of the Kondo regime, for which |e|,|e + C/|>r. 

The direct observation of the "Coulomb blockade" 
peaks and the Kondo resonance allows, unlike in previous 
studies of molecular adsorbates, for a univocal identifica- 
tion of the MO associated to the unpaired spin. As the 
spin participating to the many-body Kondo state origi- 
nates from a single orbital, we expect a well-defined value 
of Tji over the whole molecule, as confirmed by our ob- 
servations. This is in contrast with the variations of Tk 
found for Co porphyrins adsorbed on Cu(lll), which sug- 
gest a spatially-dependent, multiple orbital origin of the 
Kondo interaction, possibly due to strong intermixing of 
TM and hgand orbitals near Ep (Ref.HUI). 



C. Inelastic triplet-singlet and Kondo excitations 

In contrast to the spectra of the TM ions, the spectrum 
of NiPc recorded on benzene does not present signatures 
of finite-bias excitations. The CuPc benzene spectrum, 
however, presents two side-peaks at ±21 ± 1 meV, shown 
in Fig.[9]jb). These side-peaks are assigned to intramolec- 
ular magnetic excitations,- • owing to the presence of two 
spins in CuPc, one localized on the Cu dx2_y2 statejI^S 
and the other in the 2eg ligand state. This interpreta- 
tion is supported by the gas-phase calculation of charged 
molecules presented later. The energy of the side peaks 
provides a direct estimate of the intramolecular exchange 
coupling (J) between d- and 7r-spins. The observation of 
an intense Kondo peak at zero bias, on the other hand, 
indicates that the CuPc ground state is magnetic, namely 
a triplet (S=l) state where the metal and ligand spins are 
aligned parallel to each other, whereas the excited state 
is a singlet (S=0) state. 

The similarity between the elastic Kondo resonances 
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FIG. 9: (Color online). Fit of the differential conductance 
(dl /dV) spectra of NiPc and CuPc measured at (a) the po- 
sition of the TM ion and (b) the benzene ring, using Fano 
and Lorentzian functions for the zero and finite bias peaks, 
respectively, and step functions for the inelastic conductance 
steps. The bottom graph in (a) shows the second derivative 
(d? I / dV'^) of the spectra measured on Ni and Cu. The en- 
ergies obtained from the fits are indicated by dashed lines in 
the d?I /dV^ spectra. The bars correspond to the intensity 
of calculated Raman (yellow) and infrared (grey) modes of 
gas-phase CuPc.'^ The bottom graph in (b) shows the theo- 
retical PDOS obtained for the 2eg orbital of CuPc for 7 = 25 
meV (solid line) and meV (dashed line) calculated within 
the non-crossing approximation (see text for details). 



term (see Ref. [55] for more details) . DFT calculations 
show that the interaction between the ion big state and 
the substrate can be neglected, as reflected by the small 
hybridization-induced broadening obtained for this or- 
bital, which is a factor of 20 smaller than that of the 2eg 
state. The Hamiltonian includes a molecular exchange 
term —JSl ■ Stm, which couples the ligand spin (Sl) to 
the TM spin (Stm)- The 2eg on-site energy is fitted to e 
= -0.35 eV, and the spectral density of the coupling with 
the substrate is assumed to have rectangular shape with 
magnitude F = 0.1 eV and band half- width of 10 eV. 

The PDOS onto the 2eg state, calculated within the 
non-crossing approximation at T=5 K (see Ref. |85|, 
is compared with the experimental dl/dV spectra in 
Fig. loFb). For STM studies this comparison is a first- 
order approximation that takes into account that the 
tunneling conductance is proportional to the local DOS. 
Moreover, our calculations are performed in equilibrium. 
At typical STM parameters, the non-equilibrium modi- 
fication of NCA equationgSS for slowly varying substrate 
DOS, reduce to corrections smaller than 1/1000*'' in the 
tunneling current. Hence, in these STM studies, the equi- 
librium spectral function is an excellent approximation 
to compare with the experiment. This is confirmed by 
the correct reproduction of the Kondo features observed 
for NiPc and CuPc when we set the intramolecular cou- 
pling term to meV and 25 meV respectively, as shown 
in the bottom graph of Fig. loFb) . In the case of ferro- 
magnetic coupling, the finite-bias features correspond to 
inelastic Kondo replicas that originate from spin excita- 
tions. Decomposition into spin channels yields that the 
Stokes satellite (positive bias) is in the S = channel, 
whereas the anti-Stokes peak (negative bias) as well as 
the Kondo peak are in the S = 1 channel .I^^I^ 



of NiPc and CuPc, with S=l/2 and S=:l respectively, 
can be explained by the fact that the total CuPc spin 
is underscreened due to the very weak coupling between 
the big {dx-2_yi) orbital and the substrate. Thus, only the 
2eg ligand spin is effectively screened by the substrate in 
the temperature range accessible by our STM. 



D. Multi-orbital non-crossing approximation of the 
Kondo resonance associated to the ligand spin 

Although spin excitations coupled to the Kondo ef- 
fect have already been ob served in transport experiments 
with molecular junctionSjI^SSIllHlSl ^^j^g precise configu- 
ration of the molecule-metal bonding provided by STM 
and the assignment of the Kondo-screened spin to a spe- 
cific orbital allow us to model our results by mapping the 
electronic structure of CuPc/Ag(100) calculated by DFT 
onto an impurity Hamiltonian. Our model describes the 
interaction between the 2eg orbital and the substrate by 
an SU(4) Anderson Hamiltonian, where two-fold spin de- 
generacy coexists with the two-fold orbital degeneracy 
on the same footing, and an infinite repulsive Coulomb 



E. Magnetic moments calculated by DFT 

By calculating the saturation magnetic moments of 
TMPc, DFT provides complementary results to the anal- 
ysis of the Kondo spectra reported above. Moreover, 
in addition to supporting the interpretation of the STS 
data, DFT allows us to construct a quantitative picture 
of the effect of the different charge transfer mechanisms 
on the molecular magnetic properties. Table IH] sum- 
marizes the magnetic moments calculated for the neu- 
tral and anionic form of gas-phase TMPc as well as for 
TMPc adsorbed on Ag(lOO). We notice that, ahhough 
the amount of charge transferred from the substrate to 
the molecules is similar in all cases, the total magnetic 
moment is reduced in FePc and CoPc, whereas it is in- 
creased in NiPc and CuPc with respect to the neutral 
gas-phase TMPCs. The calculations confirm that the 
magnetic moment of the TM ions in NiPc and CuPc is 
not perturbed upon adsorption, owing to the small hy- 
bridization of the planar big orbital with the substrate 
states. This is in agreement with previous results ob- 
tained by XMCD on CuPc/Ag(100) (Ref. 32) and with 
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the coexistence of TM and ligand spins in CuPc inferred 
from the Kondo spectra. We recall that, due to the un- 
derestimation of correlation effects in DFT, the ligand 
spin is just barely present in the calculations with adsor- 
bates. For NiPc, we calculate a small magnetic moment 
of 0.14/iB, while for CuPc the moment is just slightly 
larger than the Ifig corresponding to the unpaired spin 
in the big orbital. Since charge transfer is close to one 
electron, one way to estimate the magnitude of the lig- 
and spin using DFT is to consider gas-phase anions as a 
proxy of adsorbed TMPc. In such a case the number of 
electrons belonging to the molecule is fixed and the calcu- 
lations can be constrained to yield the minimum energy 
spin state. Figure [lO] displays the spin density of neutral 
and anionic gas-phase molecules, where we clearly ob- 
serve the additional ligand spin of NiPc and CuPc upon 
charge transfer. In CuPc, the spin density originating 
from the big orbital is also observed, distributed over the 
Cu ion and Np atoms. 



responsible for the quenching of the magnetic moment. 
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1.13 
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TABLE II: Computed charge transfer AA'^ (in electrons) 
for the adsorbed molecules, and total magnetic moment m 
(units of jj,b) in neutral and anionic gas-phase, and adsorbed 
molecules. 
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FIG. 10: (Color online). Spin density of neutral (top) and 
anionic (bottom) gas-phase TMPc in D4h symmetry. Red 
(yellow) indicate spin up (down). The isovalue for the spin 
density is 0.005 e"/A^. 

FePc and CoPc present a more complicate picture com- 
pared to the simple addition of one spin in NiPc and 
CuPc, partly because of the strong hybridization with 
the substrate discussed in Section [VB| In both cases, the 
magnetic moment is reduced upon adsorption. The ex- 
tra electron donated by the surface is not easily assigned 
to a single MO, inducing a reorganization of the charge 
within the molecule. We use the results obtained for 
gas-phase anions to disentangle the different mechanisms 



As shown in Fig. II the spin density of FePc follows the 
same distribution in the gas-phase anion and in the ad- 
sorbed molecule. This similarity can also be checked by 
comparing their element-resolved magnetic moments in 
Table |III[ which are very similar in both cases. Thus, 
regardless of the strong hybridization of adsorbed FePc, 
the main features of the spin moment distribution are 
again captured by the simple addition of one electron to 
the gas-phase system. The true nature of the ground 
state of FePc on Ag (100), however, might be more com- 
plex than suggested by the calculated spin density distri- 
bution, due to the interplay of several MO close to Ep 
and dynamic electron correlation effects that are not in- 
cluded in DFT. The computed electronic structure for 
FePc yields a S=I/2 system when adsorbed on Ag (100), 
same as for the gas-phase anion. However, the adsorbate 
exhibits a more complex configuration due to the pres- 
ence of multiple molecular orbitals at the Fermi level. 
This could result in either a mixed-valence system, given 
the Fermi level crossing and strong interaction of the aig 
orbital with the substrate, or a more complex Kondo be- 
havior with a low-temperature Kondo phase originating 
in the many less-coupled orbitals near the Fermi energy. 
This complex behavior makes FePc/Ag(I00) an inter- 
esting candidate to study very low temperature Kondo 
physics in a multi-orbital system. 

CoPc represents another case of nontrivial charge 
transfer occurring upon adsorption. For example, the aig 
orbital is filled, but new empty d^ states appear above 
Ep, as shown in Figs. Kb) and (c). Moreover, the an- 
ion is not a good model for the adsorbed molecule. As 
shown in Table [ll] and Fig. [lO) the anion is in a S=0 
state, due to antiferromagnetic coupling between spins 
residing in different MOs, while in the adsorbed case 
we find a non-integer spin magnetic moment of 0.63^b. 
This non-integer spin, together with the absence of a 
Kondo interaction, could indicate that CoPc is in the 
mixed-valence regime, with charge fluctuating between 
the CoPc d and Ag states, similar to results obtained for 
CoPc on Au(III) using x-ray absorption spectroscopy.Sl 





FIG. 11: (Color online). Spin density distribution of FePc 
(a) adsorbed on Ag(lOO) and (b) in the gas-phase anion. 
Red/yellow indicate spin up/down. The isovalue for the spin 
density is 0.005 e~ / k? in both cases. 
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TABLE III: Computed element-resolved magnetic moment m 
(units of hb) for FePc in the neutral and anionic gas-phase, 
and in the adsorbate, using GGA-fU. 



F. Intramolecular exchange coupling 

The TM-ligand spin coupling observed via the inelas- 
tic Kondo interaction in CuPc can be studied in gas- 
phase [CuPc]'^". Here, the energy of ferromagnetic and 
antiferromagnetic spin configurations can be computed 
by fixing the total spin of the molecule. The ferro- 
magnetic (triplet) alignment is favored over the antifer- 
romagnetic (singlet) one. The triplet ground state is 
in agreement with the zero-bias Kondo resonance ob- 
served experimentally. From the energy difference be- 
tween triplet and singlet states, evaluated using the par- 
allel spin (E^^) and antiparallel spin (E^^) energies, and 
taking the correct singlet and triplet configurations to 
yield J = 2{E^^~E^^), we extract the exchange coupling 
constant J for the anion, J — 38 meV. The discrepancy 
with J = 21 meV measured for CuPc/Ag(100) can be 
explained by surface screening, which is expected to sig- 
nificantly reduce exchange correlation effects involving 
hybridized orbit als. 

The distribution of the spin density in [FePc]^~ and 
[CoPc]^~, shown in Fig. 10 reveals a different type of in- 



tramolecular spin correlation. Both neutral and anionic 
molecules exhibit antiferromagnetically aligned TM and 
ligand spins, indicating that the spin contrast is not di- 
rectly due to the interaction with the substrate, but is 
intrinsic to the molecule and related to the pristine TM 
spin, which extends to nearby C and N atoms upon the 
formation of MOs with mixed d- and 7r-character. Such 
intramolecular antiferromagnetic spin coupling is a result 
of the large exchange splitting of the d-levels of the Fe and 
Co ions. This induces spin-dependent mixing of the d and 
ligand states, leading to the formation of spin-polarized 
MOs with different spatial distribution. A similar mech- 
anism can be invoked to explain the site- and energy- 
dependent spin contrast recently measured for CoPc de- 
posited on Fe(llO), for which the total molecular moment 
is zero .'221 



VII. CONCLUSIONS 

The adsorption, electronic, and magnetic structure of 
four TMPC (M = Fe, Co, Ni, and Cu) adsorbed on 
Ag(lOO) have been systematically studied using STM, 
STS, and ab-initio DFT calculations within the LDA, 
CCA, GGA-fvdW, and GGA-I-U approximations. 

We find that all molecules adsorb with the macrocy- 



cle plane parallel to the surface and the TM-Np ligand 
axis rotated by ±30° with respect to the high-symmetry 
surface lattice vectors. The comparison between DFT re- 
sults using different approximations evidences differences 
in the TM-Ag interaction, but the leading role of vdW 
forces leads to similar molecule-substrate distances and 
azimuthal rotation angles. The rotation between molec- 
ular and surface symmetry axes imprints chirality to the 
frontier electronic orbitals without perturbing the struc- 
tural conformation of the molecules. This can only be 
observed in NiPc and CuPc at negative voltages, due to 
the contribution of the ai„ and partially occupied 2eg 
orbitals to the tunneling current and minor infiuence of 
the d-channels in the same energy region. 

The electronic structure of adsorbed TMPc is modi- 
fied mainly due to charge transfer from the substrate. 
Although all molecules receive approximately one elec- 
tron, charge is transferred to the ligand 2eg orbital in 
NiPc and CuPc, and to multiple MOs in FePc and CoPc, 
inducing internal charge reorganization. This difference 
is related to both energy level alignment and Coulomb 
repulsion. In addition, whereas FePc and CoPc exhibit 
strongly hybridized d_L-states near Ep, those are shifted 
down in energy in NiPc and CuPc, which display un- 
perturbed (i|| -states in this energy range. The electronic 
configuration of the Ni and Cu ions is not significantly 
affected by adsorption. 

Both charge transfer and hybridization tend to reduce 
the magnetic moment of FePc and CoPc. The computed 
electronic structure of FePc indicates a S=l/2 system 
when adsorbed on Ag(lOO), with a strong presence of 
MOs near the Fermi energy that could indicate either a 
mixed- valence system, given the strong interaction of the 
aig orbital with the substrate, or a more complex Kondo 
behavior with a low-temperature Kondo phase originat- 
ing in the many less-coupled orbitals near the Fermi en- 
ergy. The non-integer moment obtained for CoPc and the 
existence of a substantial molecular DOS at the Fermi en- 
ergy (see Fig. 4|, together with the absence of a Kondo 
resonance in this molecule, suggest a mixed-valence con- 
figuration. Opposite to the moment reduction in FePc 
and CoPc, the spin multiplicity is actually enhanced by 
charge transfer in NiPc and CuPc, due to the induc- 
tion of a ligand spin in the 2eg orbital. This transforms 
NiPc into a paramagnetic molecule and induces a triplet 
ground state in CuPc, where the ligand spin is exchange- 
coupled to the TM magnetic moment. 

The interaction between the 2eg ligand spin of NiPc 
and CuPc and the substrate gives rise to a Kondo in- 
teraction, which induces a prominent zero bias resonance 
delocalized over the perimeter of the molecules. Coherent 
coupling between Kondo spin flip and vibrational exci- 
tations induces inelastic Kondo resonances in both NiPc 
and CuPc. In CuPc, the total magnetic moment is under- 
screened as the Kondo energy scale of the TM spin is or- 
ders of magnitude smaller than the one of the 2eg ligand 
spin. However, the additional magnetic degree of free- 
dom due to the TM spin leads to the coupling of Kondo 



13 



and inelastic triplet-singlet excitations. The coexistence 
of different nonequilibrium Kondo processes related to 
vibrational and spin transitions in CuPc opens up the 
possibility to study the timescale and spatial localization 
of multiple spin excitation and relaxation channels within 
a single molecule, as reported in more detail in Ref. [71 

In general, we show that charge transfer and hybridiza- 
tion of MOs orbitals have profound effects on the elec- 
tronic configuration, magnetic moment, and transport 
properties of metal-organic complexes adsorbed on a 
metallic substrate. Such effects in TMPc do not depend 
monotonically on the electronic configuration of the TM 
ions. Rather, the energy position, symmetry, and spin 
polarization of the pristine MOs has to be considered 
within a comprehensive picture of the charge transfer 
process. Importantly, the magnetic moment of adsorbed 
TMPc can be either reduced or enhanced depending on 
the relative energy of the d± and 7r-levels, and their de- 
gree of hybridization with the substrate. Further, the 



possibility of inducing additional ligand spins upon ad- 
sorption may be relevant to establish magnetic coupling 
in extended molecular structures assembled on metal sur- 
faces. 
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